Ferromagnetic -spin glass transition induced by pressure in Gd2Mo207 
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R2M02O7 compounds show a ferromagnetic metal-insulator spin glass transition tuned by the 
radius of the rare earth ion R'^"'". We have studied Gd2Mo2 07 located on the verge of the transition, 
by neutron diffraction on a ^^"Gd isotopic sample, /iSR and X ray diffraction using the synchrotron 
radiation. All measurements were done both at ambient and under applied pressure. At ambient 
pressure, a ferromagnetic state is observed below the Curie temperature (Tc = 70 K). The ordered 
magnetic moments at 1.7 K are parallel and equal to 5.7(5) /xb and 0.8(2) /ib for Gd and Mo, 
respectively. The relaxation rate measured by /xSR evidences strong spin fluctuations below Tc and 
down to the lowest temperature (6.6 K). A spin reorientation occurs in the range 20 K<T<Tc. 
The ferromagnetic state is strongly unstable under pressure. Tc sharply decreases (down to 38 K 
at 1.3 GPa), and Bragg peaks start to coexist with mesoscopic ferromagnetic correlations. The 
ordered moments decrease under pressure. At 2.7 GPa long range magnetic order completely breaks 
down. In this spin glass state, Gd-Gd spin correlations remain ferromagnetic with a correlation 
length limited to the fourth neighbor, and Gd-Mo spin correlations turn to antiferromagnetic. The 
unique combination of three microscopic probes under pressure provides a detailed description of 
the magnetic transition, crucial for further theories. 



PACS numbers: 71.30.+h, 71.27.+a, 75.25.+Z 

I. INTRODUCTION 

Pyrochlores compounds R2T2O7, where R'^"'" and 
T^+are rare earth and transition or sp metal ions respec- 
tively, show geometrical frustration of the first neighbor 
interactions. This occurs not only for antiferromagnetic 
interactions between first neighbor Heisenbcrg moments, 
but also for ferromagnetic interactions if the neighbor- 
ing moments are constrained to lie along their local Ising 
anisotropy axes. This peculiarity of the pyrochlore lat- 
tice leads to exotic types of short range magnetic orders, 
such as spin liquids, spin ices or chemically ordered spin 
glasses, which are intensively studiedi*^. 

Among the pyrochlores, rare earth molybdates 
R2M02O7 have attracted special attention since the dis- 
covery of a crossover transition from an insulating spin 
glass state to a metallic ferromagnetic state, which can 
be tuned by the rare earth ionic radius^. The variation 
of the transition temperature which shows a universal 
curve for all studied compounds with mixed rare earth 
ions, suggests that the dominant mechanism comes from 
a change in the sign of the Mo-Mo interactions. This 
change is connected with a change in the band struc- 
ture, due to the specific energy of Mo t2g orbitals, situ- 
ated nearby the Fermi level and well separated from the 
other bands- . The insulating state observed at low rare 
earth ionic radius is attributed to the opening of a Mott- 
Hubbard gap at the Fermi level, due to strong electronic 
intrasite interactional. With increasing the rare earth 



ionic radius, Mo-Mo first neighbor interactions seem to 
evolve from antiferromagnetic type, dominated by a su- 
perexchange mechanism and frustrated by the geometry, 
to a ferromagnetic type, due to a double exchange mech- 
anism. 

In the ferromagnetic region, the ferromagnetic align- 
ment of the neighboring Mo moments favors the elec- 
tronic kinetic energy and a metallic conductivity, with 
a mechanism similar to that observed in the mangan- 
ites. The rare earth crystal field anisotropy also plays 
a role, since it yields a possible source of frustration in 
the ferromagnetic region. Such frustration seems to have 
important consequences on the conductivity properties. 
Especially, the compound Nd2Mo2 07 shows a giant ab- 
normal Hall effect at low temperatureSi^ which cannot 
be explained by the spin-orbit coupling. It is generally 
admitted that this effect is induced by the spin ice frus- 
tration of the 4f Nd^+ spins, transferred to itinerant Mo 
electrons through f-d interaction. 

Gd2Mo207 is especially interesting since the Gd ionic 
radius is situated just above the threshold for the 
metal insulator transition. This compound was thor- 
oughly investigated by many techniques. Magnetic 
measurements^-'l-*^ suggest a ferromagnetic-like transition 
at Tc, with magnetic irreversibilities occuring below. 
The Tc values reported by various groups range between 
40 - 70 K. Mossbauer measurements^" show a hyperfine 
field below the transition, and the hyperfine populations 
are found to be out of thermal equilibrium at 27 mK, 
indicating that Gd and Mo spin fluctuations are present 
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in the magnetic phase and persist down very low temper- 
ature. Heat capacity measurements^^ show two step-hke 
anomahes, suggesting two magnetic transitions, one at 
Tc (70 K) and the other well below (11 K). The transport 
properties strongly depend on sample preparation. First 
measurements on powdered samplesii show a metallic 
conductivity, with anomalies at Tq and in the 10 K - 
20 K temperature range. However, more recent data on 
high purity single crystalsi^ show an insulating ground 
state, which is very sensitive to impurity doping. In con- 
trast with Nd2Mo207, the Hall conductivity- shows the 
behavior of a conventional ferromagnet, as expected for 
Heisenberg like Gd'^^ moments (4f ) without orbital an- 
gular momentum. An investigation by Infra Red (IR) 
spectroscopy^^ shows the presence of a very small Mott- 
Hubbard gap of 20 meV, about 10 times smaller than in 
the neighboring insulating compounds (170 and 250 nieV 
for Dy and Ho, respectively). The concomitant changes 
in the optical, magnetic and transport properties suggest 
a quantum phase transition around Gd2Mo2 07, opening 
the possibility to tune the transition by pressure and/or 
magnetic field. Magnetic measurements under pressur&l^ 
and recent resistivity measurements^^ demonstrated this 
possibility at a macroscopic level. 

Up to now, there was no characterization of the type of 
magnetic order and magnetic fluctuations in Gd2Mo207. 
This is especially due to the huge absorption of natu- 
ral Gd, which makes neutron scattering experiments ex- 
tremely difficult. We present here an investigation of the 
magnetic order and magnetic fluctuations in Gd2Mo207, 
by combining neutron diffraction on a ^^°Gd isotopic 
sample with fj,SR experiments. We show that a ferro- 
magnetic collinear order is indeed stabilized well below 
Tc, but that it coexists with low temperature fluctua- 
tions. At intermediate temperatures (20 K<T<Tc) a 
reorientation of the magnetic moments occurs. 
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FIG. 1: Magnetic susceptibility measured in a static field of 
50 G, in the zero field cooling (ZFC) and field cooling (FC) 
processes. 



in Gd2Mo207 to applied pressure by neutron diffraction 
and /iSR. The pressure induced changes of the crystal 
structure and lattice constant were checked by X ray 
diffraction using the synchrotron radiation. We observe a 
pressure induced magnetic transition, which can be fully 
characterized using this unique combination of three mi- 
croscopic probes. We show that the ferromagnetic state is 
highly unstable under pressure : a rather small pressure 
of 1.3 GPa yields a decrease of the Curie temperature by 
a factor 2. At 2.7 GPa, long range ferromagnetic order 
has fully disappeared. Our results demonstrate that the 
changes in magnetism induced by pressure are equivalent 
to those induced by chemical pressure (R substitution). 
They confirm at a microscopic level the conclusions in- 
ferred from the magnetization^^. Therefore, keeping the 
same sample, a quantitative analysis of the neutron and 
muon data allows us to follow the evolution of the spin 
correlations and fluctuations throughout the transition. 

The paper is organized as follows. In section II, we de- 
scribe the sample characterization, bulk magnetic prop- 
erties and crystal structure. In section HI, we describe 
the magnetic state at ambient pressure, as studied from 
neutron diffraction and ^SR versus temperature. In sec- 
tion IV, we show the evolution of the magnetic state un- 
der pressure from X ray, neutron and /xSR experiments. 
In section V, we discuss the results in comparison with 
other experimental data and current theories. 
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FIG. 2: Neutron diffraction pattern of Gd2Mo207 measured 
at 300 K (paramagnetic phase) on the D2B powder diffrac- 
tometer, with an incident neutron wavelength of 1.594 A. Up- 
per and lower tick marks show the Bragg peak positions of the 
sample and vanadium sample holder, respectively. Solid lines 
show the calculated pattern {Rb~ 8.9%) and the difference 
spectrum. 



We also studied the sensitivity of the magnetic state 
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II. CRYSTAL STRUCTURE AND BULK 
MAGNETIC PROPERTIES 



Powdered Gd2Mo207 samples were synthesized follow- 
ing the procedure given in ref|3- We used Gd203 and 
M0O2 oxides as starting materials. The synthesis was 
made in Argon atmosphere without any excess of Mo, 
but with Ti/Zr chips to absorb oxygen traces. Samples 
were heated up to 1370 °C during 6 hours. Two anneal- 
ings at least were necessary to obtain the sample in pure 
form. A sample of 1 g with natural Gd was prepared for 
synchrotron and /iSR experiments and a sample of 0.5 
g with isotopically enriched ^^"^Gd was used for the neu- 
tron experiments. Both samples were characterized by X 
ray diffraction and magnetization. They are single phase 
and show the same magnetic transition. 

The DC magnetization (Fig. 1) was recorded in a 
SQUID magnetometer with a static field of 50 G. It shows 
a rather broad Curie transition, as already reported in 
the literature. The Tc value of 70 K reported below cor- 
responds to the inflexion point of the susceptibility curve, 
whereas the onset of the transition is situated around 80 
K. Irreversibilities of the susceptibility depending on the 
cooling process, -zero field cooling (ZFC) or field cool- 
ing (FC)-, are observed below the transition. The FC 
susceptibility linearly increases with decreasing temper- 
ature, whereas the ZFC susceptibility flattens and starts 
to drop below about 14 K. This drop suggests a blocking 
of domains walls mobility, possibly related to the spin 
reorientation discussed below. 

The crystal structure of Gd2Mo207 was investigated at 
300 K by combining powder X-ray and neutron diffrac- 
tion, the neutron pattern being measured in the high 
resolution-high flux diffractometer D2B of the Institute 
Laiie Langevin with an incident wavelength of 1.594 A in 
the high intensity version. In order to decrease the resid- 
ual absorption, the sample was placed in a hollow vana- 
dium cylindrical container. In the refinement, we made 
a specific absorption correction, taking into account the 
container geometrjii^, the neutron wavelength, and the 
isotopic composition of the sample. The linear absorption 
coefficients were estimated to = 6.3 and 4.1 cm^^ for 
D2B and D20, respectively. 

The Rietveld refinement of the D2B diffraction pat- 
tern (Fig. 2) was made using the crystallographic pro- 
grams of the FuUprof suite^^. We took the diffraction 
pattern from the vanadium container into account as a 
second phase. The refinement {Rb—8-9 %, Rf—5.9 % ) 
confirmed the structural model for a stoechiometric py- 
rochlore with space group FdSm, yielding a lattice con- 
stant a= 10.3481(2) A and an oxygen position parameter 
u= 0.3342(2) at room temperature. These results agree 
with previous determinations from synchrotron radiation 
X ray powder diffractiorii. 
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FIG. 3: Magnetic neutron diffraction pattern of Gd2Mo2 07 
measured at 1.7 K on the D20 powder diffractometer, with an 
incident neutron wavelength of 2.419 A. A pattern measured 
at 90 K ( above Tc) was subtracted. Tick marks show the 
Bragg peak positions. Solid lines show the pattern calculated 
with the coUinear ferromagnetic model {Rb= 11.3 %) and the 
difference spectrum. 



III. MAGNETIC STATE AT AMBIENT 
PRESSURE 

A. Magnetic neutron diffraction 

Neutron diffraction patterns were recorded between 1.7 
K and 90 K on the high flux diffractometer D20 of the In- 
stitute Laiie Langevin with an incident wavelength 2.419 
A, with the same sample and sample holder as for D2B. 
Absorption corrections were made as above. Magnetic 
diffraction patterns were obtained by subtracting a spec- 
trum at 90 K just above the transition. The magnetic 
pattern at 1.7 K is shown in Fig. 3. The magnetic peaks 
have hkl indices of the face centered cubic lattice. The 
(200) and (220) peaks, where there is no chemical con- 
tribution from Gd and Mo ions in the pyrochlore struc- 
ture (due to extinction from the Fd3m space group and 
special Gd-Mo positions, respectively) are also absent in 
the magnetic pattern. We refined the magnetic patterns 
using the FuUprof suitei^. At 1.7 K, we obtain a good 
refinement assuming a collinear ferromagnetic structure 
(magnetic factor Rb= 11.3 %, Fig. 3). The Mo and Gd 
moments orient along the same direction, with a ferro- 
magnetic coupling. According to the ^^^Gd Mossbauer 
measurements-^, the angle between the Gd moment and 
the <111> axis is close to 54°. This strongly suggests 
that the Gd and Mo moments lie along a <100> axis. 

With increasing temperature the intensities of the 
magnetic Bragg peaks start to decrease, but the (111) 
peak decreases much more rapidly than the high angle 
peaks (Fig. 4). It disappears around 50 K whereas the 
other peaks persist up to 70 K. Concomitantly, the mag- 
netic intensities calculated within the collinear ferromag- 
netic model start to disagree with the experimental data 
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(inset Fig. 4, Rb= 39 % at 40 K). 

The refinements exclude a zero ordered moment on 
either Gd or Mo sublattice. A global change of orienta- 
tion of both Gd and Mo sublattices is also excluded since 
it would not change the magnetic diffraction pattern of 
a powdered cubic structure. The observed anomaly is 
therefore attributed to a decoupling of Gd^+ and Mo^+ 
moments above 20 K, yielding non coUinear moments in 
the Gd (and possibly Mo) sublattice. Unfortunately, a 
precise determination of this non-coUinearity is hampered 
by the high symmetry of the crystal structure, the pres- 
ence of eight independent moments in the unit cell, and 
the dominant ferromagnetic character. 

To search for non coUinear structures, we first per- 
formed a symmetry analysis^^'"^^ using the program 
BasIrepa^Si, searching for all k=0 structures correspond- 
ing to irreducible representation of the I4i/amd space 
group, the highest symmetry group allowing ferromag- 
netism. We also checked non coUinear structures pre- 
dicted by theory, assuming either anisotropic or dipolar 
interaction322i24. Finally, we explored small deviations 
from the ferromagnetic coUinear case using the simulated 
annealing process available in FuUprof^a. Whatever the 
procedure, we could not improve the quality of the re- 
finement significantly with respect to the coUinear ferro- 
magnetic model. We conclude that a non coUinearity of 
the Gd/Mo ordered moments indeed occurs above 20 K, 
although we cannot determine its nature presently. We 
also notice a slight change of the peak linewidth, show- 
ing that the magnetic correlation length increases with 
decreasing T in this temperature range (inset Fig. 5). 
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FIG. 4: Integrated intensity of several magnetic Bragg peaks 
in Gd2Mo207, measured versus temperature on D20. The 
intensities are scaled to the intensity of the (222) nuclear peak 
at 90 K. In inset the magnetic diffraction pattern at 40 K. 

The temperature dependence of the ordered magnetic 
moments deduced from refinements in the coUinear fer- 
romagnetic model is plotted in Fig. 5. At 1.7 K, the 



moments values deduced from the refinement are 5.70(8) 
and 0.82(5) /^b for Gd and Mo, respectively. Consid- 
ering the uncertainty on the scale factor between nuclear 
and magnetic scattering due to the absorption correction, 
we estimate the ordered magnetic moments in absolute 
scale with a larger error bar, namely 5.7(5) fiB and 0.8(2) 
/iB for Gd'^"'" and Mo^"*" ions, respectively. The ordered 
Q(^3+ jjjoment is reduced with respect to the free ion 
value of 7 /iB- The ordered Mo*+ moment is also much 
smaller than the value of 2 /iB expected assuming a ionic 
description (4d^ t2g state with S=l and g=2). The cal- 
culated ordered moment per formula unit of 13(1) is 
significantly smaller than the saturated value of 16.8 /is 
obtained at 2 K in a high field of 14 T^^. Since we do 
not expect any strong reduction of the free ion values 
from crystal field effects (the Gd'^^ moment has no or- 
bital contribution) , the reduced moment values should be 
associated with the low temperature fluctuations. This 
is also supported by the /iSR experiments, as discussed 
in the next section. 

In Fig. 5a, we fitted MGd(T) by a Brillouin func- 
tion B7/2, assuming that the Gd ions are submitted to 
a molecular field coming from the Mo ions. As the Mo 
spontaneous moment is practically temperature indepen- 
dent up to Tc, we assumed that the molecular field ilex 
is constant up to 80 K and vanishes above. The fit in 
Fig. 5 is obtained with Rex = 10.8 T. 

The thermal variation of the Gd moment has also been 
inferred from the -'^^^Gd Mssbauer data in Gd2Mo207 in 
assuming that the Gd moment is proportional to 
the hyperfine field. Comparison with the present neutron 
data shows a clear disagreement: the Mossbauer derived 
values fall much more rapidly than the neutron values 
as temperature increases. So we propose a new interpre- 
tation of the Mossbauer data: the measured hyperfine 
field Hht is assumed to be the sum of a contact hyperfine 
field, proportional to the Gd moment, and of a trans- 
ferred hyperfine field Htr coming from the polarization of 
the conduction band by the Mo momentsiSiii: 



Hm{T) ^ lAMcdiT) + Ht, 



(1) 



The sign of Hhf cannot be obtained by the Mossbauer 
data. The transferred field is proportional to the Mo 
moment, and it is thus temperature independent up to 
Tq. In Fig. 5b, we have represented the hyperfine field 
data from Ref. and the fit to the above formula, 
where Mcd{T^) is calculated using the Brillouin function 
as above. The fit yields a contact hyperfine constant 
A=5.8 T//iB and a transferred field Htr= - 13 T. The sign 
of the latter is opposite to that of the contact term, which 
can occur in intermetallic Gd compounds with 3d or 4d 
metalsiSill. Above 40 K, the measured hyperfine field is 
below 4 T, with a large error bar. The above law pre- 
dicts that Hhf (T) shows a minimum at 50 K (not shown) , 
which could not be observed due to the small hyperfine 
field values. According to this picture, Gd2Mo207 is an 
example where the measured ^^^Gd hyperfine field is not 
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FIG. 5: a: Temperature dependence of the ordered moments 
Mcd (circles) and Mmo (squares) in Gd2Mo207 deduced from 
refinements of the neutron data in the collinear ferromagnetic 
model. The arrow shows the free ion value of 7 fj,B- The solid 
line for Gd correspond to a fit with a Brillouin function, as 
described in the text. In inset, the magnetic correlation length 
Lc deduced from the peak widths (dots) and the magnetic 
factor Rb (squares) versus temperature. Dashed lines are 
guides to the eye. biThermal variation of the ^^^Gd hyperfine 
field in Gd2Mo2 07, from Ref^. The line is a fit to the sum 
of a contact hyperfine field, proportional to the Gd moment, 
and of a transferred field, proportional to the Mo moment, 
and thus temperature independent. 

proportional to the spontaneous moment. 



B. /iSR measurements 

^SK measurements at ambient pressure were per- 
formed at the Swiss Muon Source at the Paul Scherrer 
Institute (Villigen, Switzerland) on the GPS instrument 
in the temperature range 6.6 K - 300 K. Selected fiSR 
spectra are shown in Fig. 6. Above Tc, the relaxation 
function Pz(t) shows an exponential like decay: Pz (t)= 
exp{-Xzt). Below Tc, Pz(t) shows a rapidly damped 



oscillation at early times which is readily attributed to 
the presence of magnetic order. An exponential decay is 
also observed at long time scales in the ordered phase, 
showing that spin fluctuations persist down to the lowest 
temperature. 
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FIG. 6: Muon spin depolarization function Pz(t) in 
Gd2Mo2 07 measured on the GPS beam line at several tem- 
peratures (ambient pressure). Lines are fits as described in 
text. The inset shows the depolarization function in an ex- 
tended time range for T=20 K. 

Below Tc, the relaxation function Pz(t) was fitted by 
the equation: 

Pz{t) = [expi^Xzt) + 2exp{~\Tt)cos{-f^ < Bioc > t)]/'S. 

This equation is expected to hold in the magnetically 
ordered phase of a powder sample??. The first term cor- 
responds to the depolarization by spin fluctuations with 
the longitudinal relaxation rate Xz, whereas the second 
term reflects the precession of the muon spin in the aver- 
age local field <B;oc> at the muon site. 7^ is the muon 
gyromagnetic ratio. The transverse relaxation rate At 
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can have both static and dynamical character. In the 
high temperature Umit, when the rate of fluctuations of 
Gd^+ and Mo^+ moments is much larger than the cou- 
pling between the muon spin and the electronic spins, one 
gets <B;oc>=0, and Xt=Xz- Then the above equation 
reduces to a simple exponential Pz{t)= cxp{-Xzt). This 
simple function was used to fit the data above Tc- 
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FIG. 7: From /xSR measurements in Gd2Mo207. a: tem- 
perature dependence of the longitudinal relaxation rate Xz; 
b: temperature dependence of 7^<B;oc> (with 7^ the muon 
gyromagnetic ratio), At and Mcd (scaled). 

In a standard ferromagnet, the longitudinal relaxation 
rate Xz shows a cusp at the critical temperature, then de- 
creases with decreasing temperature below Tc and tends 
to zero as T ^ 0. In contrast, in Gd2Mo207, Xz keeps a 
high value down to the lowest temperature (6.6 K), and 
remains T-independent in the whole ordered phase (Fig. 
7a) almost up to Tc. This abnormal behavior reflects 
the presence of strong spin fluctuations coexisting with 
the magnetic order, which seem to persist as T ^ in 
contrast with standard spin waves. Such peculiar fluc- 
tuations have already been observed by /iSR in several 



geometrically frustrated magnets with a pyrochlore lat- 
tice. They were seen in spin liquids^SiSi like Tb2Ti207 
and Yb2Ti2 07, as well as ordered compoundsSLSSiSS like 
Gd2Sn207 and Tb2Sn207. In all cases, the longitudinal 
relaxation rate saturates in the lowest T range (to values 
around Ifis^^), as here for Gd2Mo207. 

The transverse relaxation rate At is up to two or- 
ders of magnitude larger than Xz{T)- It smoothly in- 
creases below Tc and scales with the average local field 
<B;oc> as temperature varies. This leads us to assign 
At mainly to the width of the distribution of local fields. 
We notice that both quantities also scale with the or- 
dered moment MGd(T) determined by neutron diffrac- 
tion (Fig. 7b). It suggests that the local field seen 
by the muon mostly comes from the Gd'^"'' ions with 
much larger moments, although more localized, than the 
Mo^+ ones. Our current neutron and /iSR data in the 
(Tb,La)2Mo207 series'^^ as well as previous /iSR data 
in Tb2Mo207 and Y2M02O7 spin glasses'^'' support this 
interpretation, showing that the static internal field is 
about 10 times higher in Tb2Mo207 with both Tb^+ and 
Mo'*+ moments than in Y2M02O7 where only Mo''^"'" mo- 
ments are involved. 



IV. THE PRESSURE INDUCED STATE 

A. Crystal structure under pressure : X ray 
diffraction 

X-ray diffraction under pressure using the synchrotron 
radiation was performed at room temperature on the 
ID27 beam line of ESRF, in the pressure range 0-10 GPa, 
with an incident wavelength of 0.3738 A. We used a 
diamond anvil cell and an ethanol-methanol mixture as 
pressure transmitting medium. The crystal structure re- 
mains cubic with Fd3m space group in the whole mea- 
sured pressure range. 

The evolution of the lattice constant a versus pressure 
is shown in Fig. 8. From ref. 3, one estimates the critical 
lattice constant for the ferromagnetic-spin glass transi- 
tion to flc =10.327(5) A (critical ionic radius rc=1.047 
A). Besides Gd2Mo207, we also studied two samples un- 
der pressure in comparison, namely Tb2Mo207, an insu- 
lating spin glass with ao=10.312 A, smaller than Uc, and 
(Tbo.8Lao.2)2Mo207 with ao=10.378 A greater than Oc. 
In the studied pressure range, we found the same vari- 
ation of the relative volume V/Vq (where Vo is the unit 
cell volume at ambient pressure) for all compounds. The 
equation of state was fitted to the Murnaghan equation 
V/Vo = (P*Bi/Bo +l)-i/'^i, yielding a bulk modulus 
Bo=f5I.8(9) GPa . The pressure derivative of the bulk 
modulus was fixed to a reasonable value (Bi= 4.5). 

Knowing the evolution of the lattice constant under 
pressure enables us to compare the pressure induced 
state in Gd2Mo207 with the ambient pressure state in 
R2M02O7 compounds with smaller lattice constant. As 
discussed in the last section, the different experimental 
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determinations of the critical pressure for Gd2Mo207 lie 
in the pressure range 0.6-2.4 GPa. 
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FIG. 8: Variation of the lattice constant a versus pressure at 
ambient temperature. In the inset the relative volume V/Vo 
is plotted versus pressure up to 40 GPa for three samples: 
Gd2Mo207 (circles), Tb2Mo2 07 with oq < Uc (squares) and 
(Tbo.8Lao.2)2Mo207 with oo > Qc (triangles). Solid lines are 
fits to the Murnaghan equation. 
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FIG. 9: Temperature dependence of the relaxation rate Xz, 
with the Gd2Mo207 sample in the pressure cell. Data were 
measured on the GPD instrument: ambient pressure (circles) 
and P=1.3 GPa (squares). Solid lines are guides to the eye. 



T range 3.2 K -130 K, then the cell was pressurized and 
the experiment repeated. The pressure was determined 
by measuring the superconducting transition of a lead 
wire inside the pressure cell, yielding a value of 1.3 GPa. 
Below Tc it was difficult to extract any information from 
the fiSR spectra at small times (the 2/3 term), due both 
to the large background of the pressure-cell and to the 
fast depolarization of the 2/3 term. Hence, below Tc 
we fitted only the 1/3 term with an exponential depo- 
larization function. Fig. 9 compares the temperature 
dependence of the relaxation rate at ambient and under 
1.3 GPa. The sharp peak in Xzi"^) clearly moves to- 
wards lower temperatures under pressure. It reflects a 
very strong decrease of the transition temperature- from 
70 K to 38 K, when the pressure increases from ambient 
pressure to 1.3 GPa. 



We notice that the quality of the experimental data 
was not sufficient to refine the oxygen parameter due to 
texture effects and/or non isotropic powder averaging. 
So we cannot determine the pressure dependence of the 
Mo-O-Mo bond angle. 

B. Change of the transition temperature : /iSR 

Muon spin rotation measurements under pressure were 
performed on the GPD instrument of the Paul Scherrer 
Institute, using high energy incident muons to penetrate 
the pressure cell. The sample was mounted in a piston 
cylinder cell inserted in a cryostat. The signal from the 
pressure cell (about 70 % of the total asymmetry) was 
fitted by a Kubo-Gauss function. The sample was first 
measured in the pressure cell at ambient pressure in the 



C. Change of the spin correlations: magnetic 
neutron diffraction 

Neutron diffraction measurements under pressure were 
performed on the spectrometer G6-1 of the Laboratoire 
Leon Brillouin with an incident wavelength of 4.741 A, 
in the high pressure version^^. The sample was loaded 
in a sapphire anvil cell with 40 % Al powder to decrease 
the sample absorption. Measurements were performed at 
the pressures 0.5, 1.2, 1.9 and 2.7 GPa, between 1.5 K 
and 90 K. Pressure was measured by the ruby fluores- 
cence technique. The magnetic diffraction spectra were 
obtained by subtracting a spectrum at 90 K (paramag- 
netic phase). The magnetic intensity was scaled to the 
integrated intensity of the (222) Bragg peak. This proce- 
dure allows us to compare spectra measured at different 
pressures and temperatures. 
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FIG. 10: Magnetic difFraction patterns in Gd2Mo207 under 
pressure versus the scattering vector g=47rsin^/A. a. P= 1.2 
GPa; b. P= 2.7 GPa. Spectra are corrected by a spectrum 
in the paramagnetic range (90 K) and scaled to the nuclear 
intensity ol the (222) Bragg peak. 



Fig. 10a shows magnetic difFraction patterns at 1.2 
GPa. The spectrum at 1.5 K shows that magnetic long 
range (LRO) and short range (SRO) orders coexist in 
the sample. Small Bragg peaks are clearly visible, with 
a width limited by the experimental resolution. An in- 
tense scattering is also observed at low q values. This 
small angle neutron scattering (SANS), which was ab- 
sent in the ambient pressure data, shows the onset of fer- 
romagnetic correlations with a mesoscopic length scale. 
With increasing temperature, the (111) peak disappears 
at about 20 K, whereas the (222) peak vanishes at a 
higher temperature close to 38 K, the transition temper- 
ature determined by fiSR for this pressure. It suggests 
that the spin reorientation observed in the ordered state 
at ambient pressure temperature still persists under pres- 
sure. The short range ferromagnetic correlations persist 



slightly above the transition. 

At 1.9 GPa, a very small ordered component remains 
at the Bragg positions. At 2.7 GPa, the magnetic Bragg 
peaks have completely disappeared, and only magnetic 
short range order is present (Fig. 10b). Short range 
magnetic correlations are clearly shown at 1.5 K, yield- 
ing a strongly modulated magnetic background. Here, 
the SANS signal coexists with a diffuse peak centered 
around 1.1 A~^. With increasing temperature, the dif- 
fuse peak flattens and vanishes around 20 K, whereas the 
ferromagnetic correlations persist up to 40 K. 

In Fig. 11a, we compare diffraction patterns at 1.5 K 
for several pressures. With increasing pressure, the de- 
crease of the Bragg intensity coincides with the growing 
of the diffuse magnetic peak, as the relative contributions 
of SRO and LRO vary in opposite way. The SANS sig- 
nals flattens, showing that the ferromagnetic correlation 
length decreases with increasing pressure (see inset Fig. 
11a). 

To analyze the pressure data we proceeded as fol- 
lows. The long range magnetic order was analyzed within 
the collinear ferromagnetic model, by refining the Bragg 
peaks using FuUprof. In this procedure, the diffuse scat- 
tering and SANS signal were treated as a background. 
The magnetic peaks were scaled to the nuclear Bragg 
peak (222) measured in the paramagnetic range. The 
calculated spectra are shown in Fig. 11a. The ordered 
moments are plotted versus pressure Fig. lib. The 
Qj3+ moment strongly decreases with increasing pres- 
sure. The Mo'*"'' is expected to decrease too, but being 
much smaller, its variation remains within the experi- 
mental error. 

To get information about the short range correlations 
yielding the diffuse magnetic signal at 2.7 GPa, we used 
the model proposed by Bertaut and Burlet'^^ for spin 
glasses, and applied by Greedan et al.^"^ to the pyrochlore 
system. A fit of the diffuse magnetic scattering by the 
sum of radial correlation functions was performed, giving 
information on spin-spin correlation parameters 7 (Fig. 
12). We considered correlations up to the 4th neighbour 
shell (~ 7.3 A). The correlation coefficients 7 deduced 
from this fit at 2.7 GPa are plotted in inset of Fig. 12. 
The Gd-Gd correlations are F (7i,3,4> 0), while the Gd- 
Mo are AF (72 < 0). The AF Mo-Mo correlations re- 
sponsible for the frustration in the SG state2& cannot be 
detected, their contribution being about 50 times smaller 
than the Gd-Gd ones due to the smaller Mo moment. 

We outline that in the intermediate pressure range (0 
< P < 2.7 GPa) neither the ordered ferromagnetic model 
nor the short range magnetic model can describe the 
SANS signal in a proper way. This is because the meso- 
scopic length scale of the SANS is intermediate between 
the 4th neighbor shell (7 A) and the domain size associ- 
ated with the Bragg peaks (400 A or above, taking into 
account the experimental resolution). Considering the 
relative magnitude of the Mo and Gd local moments and 
the fact that Mo-Mo and Mo-Gd correlations are AF, we 
attribute the SANS signal to ferromagnetic Gd-Gd corre- 
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lations. To get an insight about their correlation length, 
we have fitted the SANS with a Lorentzian function. The 
correlation length Lc is about 35 A at 0.5 GPa and 1.5 
K. It decreases with increasing pressure (inset Fig. lib) 
and temperature. 



V. DISCUSSION 

In this section, we discuss the main results in com- 
parison with other experimental data and existing the- 
ories. At ambient pressure, the ordered ground state in 
Gd2Mo207 is a collinear ferromagnet, in contrast with 
Nd2Mo207, where the Nd^+ free ion anisotropy induces 
a spin ice configuration of the Nd'^+ moments-'^^. 

A collinear ground state is expected in Gd2Mo207 
from the band structure?-, and from the spin only char- 
acter of Gd^+ ion. It could explain why there is no gi- 
ant anomalous Hall effect, in contrast with Nd2Mo207. 
The ordered magnetic state is however strongly abnor- 
mal, even at ambient pressure. It coexists with strong 
fluctuations, which we probe by /iSR, and which persist 
down to 27 mK according to Mossbauer dataiS. These 
fluctuations could be responsible for the strong reduction 
of the ordered moments. The Gd'^^ moment is reduced 
by 20%, and the Mo^+ moment by 60 % with respect to 
the free ion values. Spin fluctuations likely play a role in 
the spin non coUinearity above 20 K. The origin of these 
fluctuations needs to be clarifled. Low temperature col- 
lective fluctuations are widely observed in geometrically 
frustrated magnets and seem to be a key feature of ge- 
ometrical frustration^-. However, in Gd2Mo207 at am- 
bient pressure, both Gd-Gd and Mo-Mo first neighbor 
interactions are ferromagnetic, and should be frustrated 
only by the anisotropy of Gd or Mo ions. This anisotropy 
is very small, and obviously does not play a key role to 
select the ground state. Alternatively, one could specu- 
late that the abnormal fiuctuations come from the metal- 
insulating instability. Taking them into account in the 
band structure is a challenge for theory. 

The ferromagnetic state is strongly unstable under 
pressure, so that a small pressure of 0.5 GPa induces im- 
portant changes in the magnetic correlations and strongly 
decreases the Tc value. Actually in Gd2Mo2 07 our neu- 
tron data show that the magnetic order gradually evolves 
from F to SG in the pressure range of 0.5 - 2.5 GPa : in 
this region, both SG and F phases coexist and their rel- 
ative amounts change with pressure. 

Initial measurements in the substituted series^ 
(RR')2Mo207 support a unique threshold from a fer- 
romagnetic metal to an insulating spin glass state for 
a lattice constant Uc— 10.33 A. Using the equation of 
state found above, it yields a chemical critical pres- 
sure of 0.7 GPa. Actually, macroscopic measurements 
on Gd2Mo2 07 under pressure show a more complex 
situation. Magnetic measurements under pressur o^^i^^ 
suggest a transition from ferromagnetic to spin glass 
state already at 0.6 GPa. On the other hand resistiv- 



ity measurementsi^ on insulating Gd2Mo207 crystals, 
show that they become metallic at a much higher pres- 
sure of 2.4 GPa. It suggests that the ferromagnetic-spin 
glass transition could be disconnected from the metal- 
insulating one. We notice that all threshold values be- 
long the pressure range of coexistence deduced from our 
neutron data. 

In Fig. 13, we have plotted the transition tempera- 
ture Tc measured in Gd2Mo207 versus the lattice con- 
stant, by combining our magnetization, neutron, muon 
and X ray data. We also plotted the spin glass transi- 
tion temperatures Tsg corresponding to the susceptibil- 
ity anomalies under pressure^. These Tgc values are in 
the temperature range 20-25 K, in good agreement with 
the values found at ambient pressure in Tb2Mo207 and 
Y2M02O7 spin glasses with smaller lattice constanljSSiSi. 
The temperature for the onset of short range correlations 
in the spin glass state is noticeably higher (around 30-40 
K at 2.7 GPa) as it often occurs in spin glasses. The 
whole set of data provides a precise description of the 
magnetic transition in the instability region. 

By comparing the magnetic pattern of Gd2Mo207 at 
2.7 GPa with patterns of R2M02O7 compounds with 
small ionic radii (R=Tb, Y and Yb) at ambient pressure, 
we can get a general insight on the different types of spin 
correlations in the spin glass phase. For this purpose we 
measured the Y, Tb and Yb compounds at ambient pres- 
sure (Fig. 14). Our results for R=Tb and Y agree with 
previous resultsSSiS^iSI. We can also follow the evolution 
of the Mo-Mo, Gd-Mo and Gd-Gd correlations with pres- 
sure from the ferromagnetic to the spin glass phase. 

In the collinear ground state found in Gd2Mo207 
at ambient pressure, all correlations are ferromagnetic. 
With increasing pressure, Mo-Mo interactions become 
AF and frustrated by the lattice geometry yielding the 
SG phase. Neutron scattering cannot probe the Mo-Mo 
correlations in the SG phase of Gd2Mo207 since their 
contribution is much smaller than the Gd-Gd and Gd- 
Mo ones due to the smaller Mo moment. But these cor- 
relations can be directly evidenced in Y2M02O7 where 
only the Mo ions are magnetic and the lattice constant 
( a=10.21 A ) is equivalent to Gd2Mo207 at 6.7 GPa. 
In Y2M02O7, AF correlations with a scale of about 5 A 
forming a four sublattice structura^^ yield a peak in the 
diffuse scattering at q= 0.44 A"-'^. A similar structure is 
observed in Yb2Mo207 ( a= 10.146 A ) where the con- 
tribution of the Yb"^+ moments (around 1 /ie) is much 
lower than for Tb and Gd. 

Concomitantly with the change in the Mo-Mo corre- 
lations, the Gd-Mo correlations become AF. The Gd-Gd 
correlations remain ferromagnetic, but their correlation 
length decreases with increasing pressure. Comparing 
Gd2Mo207 with Tb2Mo207 suggests that ferromagnetic 
R-R correlations are a general feature of the R2M02O7 
series. The length scale of these correlations is reduced 
by pressure, which enhances the AF interactions. Our 
measurements on Tb2Mo207 under pressure, to be pub- 
lished later, show that it is also the case for Tb. 
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From our data, one can evaluate the first neighbor ex- 
change interactions J and their contribution to the en- 
ergy balance. One naturally expects J^mo-Mo > J^Gd-Mo 
> i/cd-Gd- At ambient pressure, the three interactions 
are likely all ferromagnetic. Dominant Mo-Mo interac- 
tions stabilize the ferromagnetic order but may be sen- 
sitive to temperature due to the band structure instabil- 
ity. Gd-Mo ferromagnetic interactions account for the F 
alignment of Gd and Mo moments. Gd-Gd interactions 
cannot be measured. One can speculate that their fer- 
romagnetic character stabilizes the collinear state at low 
temperature. 

We first estimate the exchange energy per ion (Mo or 
Gd) for the three interactions. For Mo-Mo exchange, one 
can state that Emo-Mo ~ Tc ~ ^cw = -70 K where 6'cw 
is the Curie- Weiss constant. For Gd-Mo exchange, we use 
the exchange field acting on Gd, determined by fitting 
the temperature dependence of the ordered Gd moment 
(Fig. 5): = 10.8 T; then: Ecd-Mo^ -^Gd^.. - 
-45 K. Finally, to estimate the Gd-Gd exchange/dipole 
energy, we use the ^cw value in Gd2Ti2 07 or Gd2Sn207, 
where there is no 3d/4d magnetic metal ion: Ecd-Cd ~ 
-10 K. Thus Gd-Mo and Gd-Gd interactions may play a 
role to determine the ground state in the ferromagnetic 
phase and threshold region. 

In order to estimate the exchange integrals, we make 
the following assumptions: the exchange between a pair 
of ions A and B is written Fox = - J a-b SaSb, where 
Sa and Sb are the true spins of the A and B ions. Then, 
the exchange energy per ion is Ea-b= -zJ^a-b SaSb, 
where z is the number of nearest neighbors (z=6 in the 
pyrochlore lattice for all exchange bonds). Although an 
ionic description is not fully appropriate for the Mo ion 
in Gd2Mo207, we assume S=l for Mo, and for Gd'^^, we 
use S—1/2 as there is no crystal field splitting. Then, we 
obtain the following exchange integrals: JT'mo-Mo'^ 12 K, 
Jod-Mo ~ 2 K and Jcd-Gd ~ 0.14 K. 

Under applied pressure, we find that Mo-Mo ex- 
change interaction changes sign and become AF. Cal- 
culated phase diagram for given values of the exchange 
interactions24j show that AF Mo-Mo interactions may 
stabilize a degenerated ground state whatever the sign 
of the Gd-Mo interactions. Mo-Mo interaction involves 
Mo t2g orbital, which are very sensitive to the electron 
correlations. Pressure is expected to increase the intra 
site electron correlations energy U. At ambient pressure, 
our estimation of J7'mo-Mo(12K) has the same order of 
magnitude as the theoretical estimation of about 5 meV 



(58K), obtained by SolovyevS- in the mean field Hartree- 
Fock approach, taking the on-site Coulomb interaction U 
=2.5 eV between Mo(4d) electrons. The model predicts 
that Jmo-Mo decreases with increasing U, and one could 
speculate that in Gd2Mo207, the J7mo-Mo exchange in- 
teraction becomes AF for a U value of about 3.5 eV, as it 
occurs in Y2M02O7 for a lower U value. Gd-Gd interac- 
tions between localized Gd moments are expected to be 
less sensitive to pressure. 

The change the Mo-Mo interaction is connected with 
the aperture of a Mott-Hubbard gap in the Mo t2c, band, 
so that the theory^ predicts an insulating spin glass 
state. The is actually true for the spin glasses with small 
ionic radii (R=Y, Ho, Dy) but not for the pressure in- 
duced spin glass state. Resistivity measurements ^^'34 on 
Gd2Mo207 and Sm2Mo207 show that the spin glass state 
induced under pressure is actually metallic. With re- 
spect to chemical pressure, an applied pressure should 
not only increase the electron correlations responsible 
for the Mott-Hubbard gap, but also increases the band 
width. This should naturally favor electron delocaliza- 
tion as a dominant feature. Optical measurements could 
also check the evolution of conduction properties under 
pressure. 



VI. CONCLUSION 

We have observed an abnormal ferromagnetic phase 
with strong fluctuations in Gd2Mo207. By applying 
pressure, we can tune the change of this ferromagnetic 
phase into a spin glass phase. The fact that an ap- 
plied pressure is equivalent to a chemical pressure to 
induce the ferromagnetic -spin glass transition supports 
a mechanism mostly controlled by Mo-Mo interactions. 
The combination of three microscopic probes under pres- 
sure allows us to follow the evolution of magnetism with 
the lattice constant with great precision throughout the 
threshold and to evaluate the role of the rare earth inter- 
actions in the energy balance. 
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FIG. 11: a. Magnetic diffraction patterns measured in 
Gd2Mo207 at 1.5 K for several pressures. Solid lines are 
fits of the Bragg peaks (0.5, 1.2 and 1.9 GPa, q > 0.5 A"^) 
and of the diffuse peaks (2.7 GPa), respectively. The dashed 
lines are fits of the SANS signal (0.5, 1.2 and 1.9 GPa, q < 
0.5 A"^). All fits are described in text. b. Variation of the 
ordered moments versus pressure. The arrows indicate the 
critical pressure as determined from susceptibility^'' (0.6 GPa, 
grey), chemical pressure- (0.7 GPa, white) and resistivity— 
(2.4 GPa, black). In inset: the variation of the correlation 
length versus pressure as obtained from the fit of the SANS. 




q (A-^) 

FIG. 12: Magnetic intensity in Gd2Mo207 at 1.5 K and P=2.7 
GPa. The fit is made with the short range model described 
in the text. In inset the temperature dependence of the cor- 
relation coefficients. Dashed lines are guides to the eye. 
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FIG. 13: Curie and spin glass transition temperatures Tc and 
TsG versus lattice constant in Gd2Mo2 07: (crosses) from sus- 
ceptibility, (squares) neutron diffraction and (diamonds) ^SR 
(this work); (triangles) from susceptibility under pressure—. 
The TsG values of Y2M02O7 (a= 10.21 A) and Tb2Mo207 
{a— 10.31 A) from references—*— are plotted for comparison. 
The arrows indicate the critical lattice constant as determined 
from susceptibilitji^ (grey), chemical pressure^ (white) and 
resistivit}"— (black), using our compressibility data. 
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FIG. 14: Magnetic intensity of R2M02O7 (R= Gd, Tb, Y 
and Yb) versus the scattering vector q. A spectrum in the 
paramagnetic phase was subtracted. 



